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Abstract 

Biopolymers serve as one-dimensional tracks on which motor proteins move to perform their biolog¬ 
ical roles. Motor protein phenomena have inspired theoretical models of one-dimensional transport, 
crowding, and jamming. Experiments studying the motion of Xklpl motors on reconstituted an¬ 
tiparallel microtubule overlaps demonstrated that motors recruited to the overlap walk toward the 
plus end of individual microtubules and frequently switch between filaments. We study a model 
of this system that couples the totally asymmetric simple exclusion process (TASEP) for motor 
motion with switches between antiparallel filaments and binding kinetics. We determine steady- 
state motor density profiles for fixed-length overlaps using exact and approximate solutions of the 
continuum differential equations and compare to kinetic Monte Carlo simulations. Overlap motor 
density profiles and motor trajectories resemble experimental measurements. The phase diagram 
of the model is similar to the single-filament case for low switching rate, while for high switching 
rate we find a new low density-high density-low density-high density phase. The overlap center 
region, far from the overlap ends, has a constant motor density as one would naively expect. How¬ 
ever, rather than following a simple binding equilibrium, the center motor density depends on 
total overlap length, motor speed, and motor switching rate. The size of the crowded boundary 
layer near the overlap ends is also dependent on the overlap length and switching rate in addition 
to the motor speed and bulk concentration. The antiparallel microtubule overlap geometry may 
offer a previously unrecognized mechanism for biological regulation of protein concentration and 
consequent activity. 
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Introduction 

The motion of motor proteins on biopolymers is important for diverse biological processes ©• 
Actin, microtubules, and nucleic acids can serve as one-dimensional tracks on which motor proteins 
(including myosins, kinesins, helicases, and ribosomes) move (O [3]). Motors must accumulate on 
filaments in sufficient density to perform their biological roles. This motor accumulation along a 
biopolymer is affected by two key effects: the directional walking of motors and motor binding 
to/unbinding from the filament. 

Motor accumulation on filaments is related to extensive theoretical work on asymmetric exclu¬ 
sion processes (ASEP) (jj]). In ASEP models, particles move on a one-dimensional (ID) lattice by 
biased hopping and experience excluded-volume interactions with other particles. These models 
have been applied to diverse examples of one-dimensional nonequilibrium transport ranging from 
molecular motors to vehicular and pedestrian traffic. Active particle motion leads to non-zero flux 
of particle density and nontrivial flux and density profiles. To compare to experimental studies of 
molecular motors, ASEP models have been extended to incorporate important biophysical ingredi¬ 
ents such as binding kinetics. A model extending the totally asymmetric exclusion process (TASEP) 
to include motor binding and unbinding (Langmuir kinetics) predicted motor density profiles along 
single fixed-length filaments ([2]). Experimental work measured kinesin-8 motor protein traffic jams 
on microtubules and found good agreement with the predicted density profiles ([6|) . 

Motor motion on cytoskeletal filaments is important for biological length regulation, including 
regulation of the length of the polymer (jTI), lengths of overlap regions between filaments (j8)), and 
the length of cytoskeletal assemblies such as the mitotic spindle daiinj) and even whole cells dm 
[ED. The case of regulation of microtubule (MT) length has seen the most work. MTs undergo 
nonequilibrium polymerization dynamics characterized by switching between distinct growing and 
shrinking states. While this dynamic instability alone leads to a broad distribution of MT lengths 
dm, numerous proteins targeted to MT ends modify their dynamics and can dramatically alter 
the length distribution (jl4|] . Single MTs can have their length regulated, for example, by kinesin-8 
motors that walk with directional bias and shorten the MT from its end dzmaiisD- Theoretical 
work on TASEP-like models has described how length-dependent depolymerization affects otherwise 
static filaments danziiisj), filaments with simplified polymerization kinetics, (jl9H22p . and dynamic 
MTs dssHm- 

Because the mitotic spindle includes arrays of overlapping antiparallel MTs at the spindle 
midzone, regulation of MT overlaps is important for mitotic spindle function and cytokinesis. The 
MT crosslinking protein PRC1/Asel/MAP65 and kinesin-4 motors (chromokinesins) play roles in 
maintenance of the spindle midzone in anaphase (I26H28P , along with other motors and MAPs ()29l) . 
Direct binding interactions of PRCl and kinesin 4 can cause length-dependent accumulation at 
plus ends of single MTs ()30p . Bieling, Telley, and Surrey (BTS) reconstituted a minimal system 
of stable antiparallel MT overlaps in which PRCl bound preferentially to overlapping regions of 
antiparallel MTs ([8|) (see also Subramanian et al. (l3ip L PRCl recruited the kinesin-4 motor Xklpl 
to the overlap. Xklpl motors could bind to and unbind from the MTs, walk toward the plus end of 
each MT, and switch between the two MTs at a relatively high rate ([8]). Motors present near the 
MT plus ends slowed the polymerization speed, consistent with earlier work showing that Xklpl 
inhibits dynamic instability (|32p and affects spindle MT mass (j33p . As a result, antiparallel MT 
overlaps reached a constant length that depended on the bulk concentration of motors. This work 
demonstrated that motor-dependent regulation of dynamics and length can occur not just for single 
MTs, but for overlapping MT pairs. 

Here we model the motor density profiles of motor proteins on antiparallel MT overlaps of fixed 
lengths. We do not explicitly consider MT length regulation in this system, as the combination 
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of binding kinetics and coupled switching of the motors between the two antiparallel filaments is 
sufficient to produce a rich phenomenology. Our work is an extension of Parmeggiani, Franosch, 
and Frey’s TASEP with binding kinetics on a single filament ([5l l3^ to include two antiparallel 
filaments coupled by switching; it is also an extension of TASEP models of two antiparallel lanes 
with lane switching ()35l l36|) to incorporate binding and unbinding kinetics. 

We first develop the model and show that measured and estimated parameters can give overlap 
motor density prohles and motor trajectories qualitatively similar to those found in the BTS ex¬ 
periments. We then develop an analytic solution to the mean-field steady state equations, and use 
it to determine the phase diagram of the model. For high motor switching rate between filaments, 
we find a new low density-high density-low density-high density phase. We then study the model 
for the reference parameter set in more detail. Because the motor density profiles are controlled 
by both boundary conditions and a non-local total binding constraint, the density in the center 
of the overlap depends not just on motor binding kinetics but also on motor speed and overlap 
length. We find an analytical approximation that describes the overlap center density and the size 
of the motor-dense boundary layer near the overlap ends. The degree of motor accumulation near 
the overlap ends depends on the overlap length and motor switching rate in addition to the motor 
speed and bulk concentration. The coupling of motor motion, binding, and switching kinetics on 
antiparallel MT overlaps may therefore offer a previously unrecognized mechanism for the control 
of one-dimensional motor density profiles. 


Materials and Methods 

Model 

In this section, we develop a mathematical model of motor density on antiparallel filament overlaps, 
inspired by the experiments of Bieling, Telley, and Surrey (BTS) (jSj). In the BTS experiments, 
the crosslinking protein PRCl binds preferentially to overlapping regions of antiparallel MTs and 
recruits the kinesin motor protein Xklpl to the overlap region. Because PRCl and therefore 
the motors are present at much higher concentrations in the overlap, in the model we consider 
the overlap region only and consider the regions of single filaments as sources or sinks of motor 
proteins (fig. 0A)- While previous work has shown that PRCl/Asel alone can develop density 
inhomogeneities and exert forces on sliding MTs (I37( I38p . in the BTS experiments little to no MT 
sliding occurred and the PRCl distribution was uniform in the overlaps ([8[). Therefore we do not 
explicitly model the PRCl molecules or their spatial distribution, but assume they are uniformly 
distributed so that each lattice site in the overlap is identical with a binding affinity that would 
correspond to the density-weighted average of the PRCl and bare tubulin affinities. We treat each 
MT as a single track and neglect the multiple-protofilament structure of the MT, consistent with 
previous theoretical work CZlEMaEl. In our model, motors can bind to and unbind from each 
of the MTs, walk toward the plus end of a MT, and switch between the two MTs (hg. Hh)- 

There are two competing processes in this model: motor stepping (TASEP) and motor binding, 
unbinding, and MT switching (Langmuir kinetics). The Langmuir kinetics dominate the system 
behavior if the overlap length is sufficiently large, N 1, so that each bound motor can only 
walk for a short fraction of the overlap length before it unbinds. Thus, the competition between 
TASEP and Langmuir kinetics occurs only if the overall binding rate to one MT in the overlap 
KonC = NkonC (where c is the bulk motor concentration), and unbinding rate ATog = Nkos, are of 
similar magnitude to the motor speed v m- 

In the next section, we develop the discrete microscopic model and corresponding stochastic 
simulation, then derive the mean-held continuum description. 
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Discrete model 

We consider two overlapping antiparallel MTs of fixed length, so the number of sites N is fixed 
(% [i^)- At each site, motor binding or unbinding can occur with binding rate konC, where /con 
is the binding rate constant per site and c the bulk motor concentration, and unbinding rate kos- 
Each bound motor steps at rate v to the next site toward the MT plus end (if the next site is 
unoccupied), and switches at rate s to the site on the adjacent MT (if that site is unoccupied). 
Note that Xklpl is plus-end directed, and we therefore assume that the motion on a single MT is 
unidirectional. 

The occupation number hi is 1 if site i is occupied or 0 if site i is empty. Then the equations 
for interior sites (2 < i < X — 1) on MTs with plus end to the right (R) and left (L) are 

= vhR^i-i{t)[l - hR^i{t)] - vhR^i{t)[l - hij,i+i(t)] + fconc[l - nR,i{t)] 

-koshR^i{t) - shR^i{t)[l - hL,i{t)] + shL,i{t)[l - hR^i{t)], (1) 

= vhL,i+i{t)[l - hL,i{t)\ - vhL,i{t)[l - hL,i-l{t)\ + konc[l - hL,i{t)] 

-koshL,i(t) - shL,i (t)[1 - UR^i (t)] -b suR^i{t)[l - hL,iit)]. (2) 

At the boundary sites, we modify these equations to incorporate fluxes into and out of the overlap. 
The flux into the overlap Is va[\ — hi(t)], where nonzero a results from motors moving into the 
overlap from the adjacent single-MT region. The flux out of the overlap is v/3hi\f{t), where /3 is 
derived fromthe rate at which motors at an MT plus end unbind. In principle, each filament could 
have different boundary conditions. Because we have no physical reason to distinguish the two 
halves of the overlap, we focus on the symmetric case or = oir = a and Pl = Pr = P- The 
resulting boundary site equations are 

va[l - hR^i{t)] - vuR^i(t)[1 - hi ?,2 (t )], (3) 

vhR^N-iit)[l - hi?,Ar(t)] - vf3hR^N{t), (4) 

vhL,2{t)[l - hL,i{t)] - vl3hL,i{t), (5) 

va[l - hL,Nit)] - vhL,Nit)[l - hL,N-l(t)]. (6) 

As shown below, these boundary conditions fix the motor densities to be a at the minus end and 
1 — /3 at the plus end of each filament. 

At steady state, we can derive a total binding constraint on the equations by summing all of 
sites on both filaments. The bulk flux terms of the form hj_i(t)[l — hj(t)] sum to zero and only the 
binding and boundary terms remain: 


dhRj(t) 

dt 

dhR^Njt) 

dt 

dhL,iit) 

dt 

dhL,N{t) 

dt 


N-l 

^ [2konC - (konC + fcofr)(hi?,i + hL,i)] + va[2 - hi?,i(t) - hL,N{t)] - vf3[hR^N{t) + nL,i{t)] = 0. (7) 
i=2 

We find a binding constraint on the total motor binding 


N 


'^nR,i + nL,i = 2N 


k„ 


2=1 


konC- feoff 


■ + 


2v[a{l -a)- /3(1 - /3)] 


feonC “h fe< 


off 


= 2Npo + 


2v[a{l-a) -/3(1 -/3)] 


feonC fe( 


off 


( 8 ) 
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where we have defined the Langmuir density po = konc/{konC + kog). Therefore, at steady state 
an equilibrium involving binding, unbinding, and the filament-end boundary conditions must be 
reached on average for the entire overlap. This is related to the zero-current condition found in 
previous work on the antiparallel TASEP without binding kinetics (j35| l36|) . We did not find an 
analytical solution to the discrete equations. Instead, we performed kinetic Monte Carlo (kMC) 
simulations of the discrete model (Supporting Material) and compared to solutions in the continuum 
limit. 

Mean-field continuum model 

We derive the mean-field continuum model as in previous work ([5]) by taking the stationary average 
{hi) = Pi, applying the random phase approximation (hihj+i) = (hj)(hi+i). We also assume 
motor commutation during track switching {nR^ihi^i) = {hi^i) {hji^i) to give discrete mean-field 
equations with linear switching terms (Supporting Material). We take the continuum limit and 
nondimensionalize the parameters and variables by choosing the length of the overlap, L, as the 
unit of length and L/v as the unit of time. Capital letters denote the nondimensionalized parameters 
(5 = sL/v and so on. Supporting Material). We choose x = 0 as the center of the overlap, and 
the boundary conditions become /?i?(— 0.5) = /9l(0.5) = a and pji{0.5) = pl{—0-5) = 1-/3. The 
steady-state continuum mean-field equations are 

0 = {2pji - + Konc{l - pr) - KosPr - SpR +SpL, (9) 

0 = {l-2pL)^+KonC{l-pL)-K,gpL + SpR-SpL. (10) 

In the continuum mean-field model, the total binding constraint (Eqn. is 

KonC a{l - a) - I3{1 - /3) _ «(! - a) - /3(1 - /3) 

KonC + KoS L{KonC + Kos) “ L{KonC + K^s) 

(11) 


ri/2 ri/2 

/ dx PR = 

1-1/2 J-1/2 


Results 

Motor density and trajectories in overlap 

To study whether our model can qualitatively describe the motor density in antiparallel MT over¬ 
laps observed experimentally, we determined a reference parameter set corresponding to the BTS 
experiments (table [^. Most parameters were directly measured by BTS. We estimated the motor 
on rate constant kon by comparison to the single-molecule imaging of low-density Xklpl in an over¬ 
lap (j8|). To estimate the minus-end boundary condition a, note that BTS found that Xklpl binding 
was greatly increased on overlaps due to recruitment by PRCl. In the BTS experiments, low ionic 
strength (which favors motor-MT binding) was needed to observe significant Xklpl binding to sin¬ 
gle MTs outside of overlaps. Therefore, we assume that the flux of motors into the overlap from 
outside is negligible, and set a = 0 in our reference parameters. The plus-end boundary condition 
/3 is related to the motor unbinding rate at MT plus ends. While this end off rate was not directly 
measured by BTS, kinesin motors typically pause at MT plus ends and have an end unbinding 
rate smaller than the unbinding rate in the bulk. We therefore expect that f3 lies between 0 (no 
end unbinding) and 2.7 x 10“^ (the value of (3 corresponding to an unbinding rate equal to kog, 
the motor unbinding rate in the bulk). We found that (3 = 2.7 x 10“^ is so small that the motor 
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Symbol 

Parameter 

Reference value 

Notes 

V 

Motor speed 

0.5 pm s”’^ 

Measured by Bieling et al. (|S]); varied up to 

8 pm s“^ to study effect of varying speed 
on density profiles 

h 

''■on 

Binding rate constant 

2.7 X 10-^ nM-i s“^ 

Estimated based on motor density profiles 
and kymographs to give affinity of 6.3 nM, 
~15x higher than measured for motor on 
single MT by Bieling et al. ll551l 

C 

Bulk motor concentration 

1-200 nM 

Varied by Bieling et al. ([81 

^off 

Unbinding rate 

0.169 s"^ 

Measured by Bieling et al. (|8| 

S 

Switching rate 

0.44 s“^ 

Measured by Bieling et al. (|H1) 

a 

Motor flux constant into 
overlap from MT minus end 

0 

Motors bind primarily inside the overlap; 
see discussion in the Results. We varied a 
between 0 and 1 to determine the model 
phase diagram 

d 

Motor flux constant out of 
overlap from MT plus end 

0 

An upper bound on the end motor unbind¬ 
ing rate is f 5 = 2.7 x 10“^; see discussion 
in the Results. We varied /3 between 0 and 

1 to determine the model phase diagram 

N 

Number of sites 

120-2500 

Varied to study overlaps of length 1-20 pm 

5 

Length of a single site 

8 nm 

Length of an a-jS tubulin dimer, see Bray 

0 


Table 1: Parameter values for the reference parameter set, taken from experimental measurements 
or estimated as noted. 


density profiles were indistinguishable from those with /3 = 0 (hg. |S1[ ) . Therefore we used /3 = 0 in 
our reference parameter set. 

With these parameters, we studied kMC simulations of our model with varying motor concen¬ 
tration and overlap lengths corresponding to the steady-state values measured by BTS. We made 
simulated images that represent how our model’s motor distributions would appear in an experi¬ 
ment with fluorescently tagged motors (green) and antiparallel MT overlap region (red) (fig. [^, 
Supporting Material). Motors decorating the overlap, with greater density and end accumula¬ 
tion for higher bulk motor concentration. The simulated images are qualitatively similar to the 
experimental images (fig. 4 in Bieling et al. (j 8 |)). We also made simulated kymographs (fig. Ep. 
Supporting Material), which show directed motor motion with direction reversal, qualitatively sim¬ 
ilar to the experimental results (hgs. 6 and S5 in Bieling et al. dS])). We note that the experimental 
kymographs show a larger population of paused/immobile motors than occurs in our model. 


Analytic solution of the steady-state continuum equations 


To further study the behavior of our model, we determined analytic steady-state solutions of the 
continuum mean-field equations. One solution to Eqn. ido] is the constant Langmuir density set 
by binding/unbinding equilibrium, po = Koncj(KonC+Kos). To hnd spatially varying solutions, we 
first dehne (Th^l = Pr,l — 5 , the difference of the motor occupancies from |. The rate combinations 
are k = KonC + Kqs + S and 7 = Ko^c — Kos- Then the equations become 


daR 

dx 

dai 


k 7 Sai 

2 Aor 2aR 

k 7 Sapt 

- ^ - - - ^ - -■ 

2 4(Jl 2aL 


( 12 ) 


dx 


( 13 ) 
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Eqns 12 and 13 are well defined for 7^ 0. This allows solution of the differential equation 
relating aji + ai and or — ul with one integration constant (Supporting Material). Our kMC 
simulation results for motor density in the overlap differ substantially from this analytic solution 
(fig. [82]), because the analytic solution (Eqn. S18) does not satisfy the total binding constraint 
(Eqn. |ll|). Therefore, the full density profiles must be determined by matching continuum solu¬ 
tions corresponding to different integration constants. This mathematical result implies several 
remarkable properties of the density profiles, as discussed below. 


Phase diagram 

We determined the phase diagram of the model as a function of the boundary condition parameters 
a and /3, as well as how the phase diagram changes with the other model parameters. We find four 
phases previously observed for the single-filament case (j5|), the low density (L), high density (H), 
low density-high density (LH), and Meissner (M) phases (figs. [2, S3). We distinguish two different 
classes of L, H, and LH phase on overlaps: the center-accumulating case (abbreviated c) has a 
higher total motor concentration profile near the center of the filament and motor accumulation 
at the overlap center, while the end-accumulating case (e) has a lower total motor concentration 
profile near the filament center and motor accumulation at the overlap ends. In addition, for high 
switching rate we observe a new low density-high density-low density-high density (LHLH) phase 

(fig- HI). 

We studied flows in the ctr — gr phase plane to determine the motor concentration profiles 
(fig. [S4) ) and therefore the phase diagram, as discussed in the Supporting Material. Representative 
phase diagrams for small and large switching rate are shown in fig. We note that the new 
LHLH phase occurs when a < 0.5 and the switching frequency is sufficiently high, conditions which 
apply to the BTS experiments. This phase can appear for high switching rate when two analytic 
solutions to the steady-state continuum equations intersect, causing transition points in the phase 


plane where the density is ill defined (fig. S4-S5). Our LHLH phase is reminiscent of multi-phase 


coexistence found by Pierobon (I40p . However, in this previous work the multi-phase co-existence 
occurs due to a point defect in the lane, while here it occurs for spatially uniform dynamics. 

The phase boundaries change with parameters of the model. Increases (decreases) to the Lang¬ 
muir density (through changes in bulk motor concentration or binding/unbinding rate constants) 
shift the lower boundary of the LH phase closer to (farther from) the line where 1 — /3 = 0.5 when 
a> pq and closer to (farther from) the line a = f3 when a < pQ. Similar shifts occur for the upper 
boundary of the LH phase. The motor speed affects the width and the shape of the LH phase: 
for higher motor speed, the LH phase region becomes narrower in width and the LH boundaries 
become more flat. 


Density profiles for reference parameters 

The motor density profiles that we determine by kMC simulation are qualitatively consistent with 
those observed experimentally (figs. 00. Here we illustrate our results for the reference parameter 
set (table ; in the Supporting Material we also discuss a large-5 parameter set chosen to approx- 

fig- 1 


imate the limit of large switching rate (table SI 


This parameter regime typically shows 
the LH phase, where the density profiles have three regions separated by domain walls. For motors 
moving to the right, there is a boundary layer on the left with the motor density increasing from 
zero, then a region of approximately linearly varying density, then a sharp transition to another 
boundary layer of linearly increasing density that approaches 1 on the right boundary (fig. [^. If 
we fix all other parameters, the boundary layer regions increase in size as the bulk motor concern 
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tration increases. The total motor density in the overlap is the sum of the motor densities on the 
two filaments, which have the symmetry pl{—x) = pr{x). The linear density variation near the 
overlap ends can be seen in the continuum mean-field equations (Supporting Material): the slope 
is KonC + 5 at the minus end and Kqs -|- 5 at the plus end. This approximation agrees well with 


simulation results near the overlap ends (fig. S2 S7). 


Control of center density by the total binding constraint 

Motor density profiles must satisfy the total binding constraint of Eqn. [H which requires that 
when a = P = 0, the integral of the density on a single MT must equal pQ, the Langmuir density 
determined by binding/unbinding equilibrium. We verified that the total binding constraint is 
satisfied in our simulations by determining the integrated motor density and comparing it to po 
(fig. [S8|). However, as noted above, the analytic solutions we found for the steady-state motor 


density (Eqns. S18, S23, S24) do not in general satisfy the total binding constraint. As a result, 


the motor concentration in the center of the filaments does not necessarily approach the Langmuir 
density po = ATonc/(ATonC + A'ofr) (fig- g. This is a significant difference between the model of 
microtubule overlaps we study and the single-filament model ([5]) . 

The total binding constraint also controls the length of the boundary layer regions near the 
overlap ends. To derive an analytic approximation expression for this length, we approximate the 
motor densities as piecewise linear, and assume the domain walls are infinitely thin so that we can 
neglect them in integrating the concentration (fig. S9). A filament is divided into three regions: 
boundary layers near the plus and minus ends, and the central region. We can then determine the 


motor density profile (Eqn. S28), the boundary layer ends, xu (Eqn. S30), and the motor density 


in the center of the overlap, pc (Eqn. S31). 


The motor density at the center of the overlap pc^ po, even when the overlap is long. Instead 
the center density depends on the motor speed and filament switching rate in addition to binding 
parameters (fig.|^ In fig. we compare simulation results and the prediction of Eqn. S31 showing 
that our approximation is in good agreement with simulations. (Note that our approximation of 
the density as linearly varying becomes less exact for slow motor speeds, leading to a discrepancy 
between simulation results and our analytic approximation.) For comparison, we show the Langmuir 
density po for the same parameter. Varying motor speed has no effect on po (right panel), but a 
significant effect on p^. To further illustrate the dependence of the center density on motor kinetic 
parameters, we show in fig. [^examples of motor density profiles from kMC simulations with varying 
motor speed. Changing motor kinetic parameters can have a large effect on the center density in 
the overlap, even for systems with identical binding kinetic parameters and po- 


Control of end accnmnlation by the total binding constraint 

The length of the boundary layer at the overlap ends depends on motor binding and kinetic pa¬ 


rameters as well as the overlap length (Eqn. S30). We compared the analytic approximation to 
simulation results and found good agreement (fig. SIO), allowing us to predict how varying experi¬ 


mental control parameters will alter the degree of motor accumulation near the overlap ends (fig. [^. 
We note that although the boundary layer lengths predicted for typical experimental parameters 
are below the optical resolution limit of ~ 250 nm, they could be distinguished by fits of fluorescent 
motor intensity profiles to model predictions (c.f. figs. [I10i , similar to how fits to a Gaussian 
fluorescence intensity profile allow sub-resolution localization of single molecules ()41l) . 

The total binding constraint means the the boundary layer length is determined as fraction of 
the overlap length. In the single-filament model, the motor density profile in the MT minus-end 
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boundary layer is independent of MT length (HZUMl), with important consequences: MT length 
regulation by the kinesin-8 motor Kip3 ([71 fTHl) depends on a motor density profile that maintains 
the same functional form near the minus end as MT length changes. By contrast, on antiparallel 
MT overlaps the density profile, boundary-layer length, and extent of motor accumulation at the 
overlap ends change with overlap length (fig. . This means that a length regulation mechanism 
identical to that of Kip3 cannot occur on antiparallel MT overlaps. 


Conclusion 


We studied a model of motor motion on antiparallel MT overlaps that incorporates motor binding 
and unbinding, plus-end directed motor motion, and switching between filaments (fig. [^. Our 
model is inspired by the experiments of Bieling, Telley, and Surrey on the motion of the kinesin- 
4 motor Xklpl on antiparallel MT overlaps (|S|). Our model is an extension of previous theory 
that studied motor motion on a single MT with binding kinetics ([5l [M|) , or motor motion on two 
antiparallel filaments with lane switching, but no binding kinetics 033 EH])- To our knowledge, this 
is the first theoretical study of a two-lane TASEP model with oppositely oriented lanes, switching, 
and binding/unbinding. 

To compare to the BTS experiments, we used measured or estimated parameters (table [^. 
Because PRCl recruits motors directly into the overlap and motor binding to single MTs is much 
weaker, we neglect motor binding to MTs outside of the overlap and set the minus-end flux param¬ 
eter a = 0. Using the bulk motor unbinding rate kos as an upper bound on the plus-end unbinding 
rate, we find /3 < 2.7 x 10“^. This value is sufficiently small that it gives motor density profiles 


indistinguishable from those with /3 = 0 (fig. SI). Therefore no-flux boundary conditions with 
a = /3 = 0 are used to model the experiments. Simulated images of the motor distribution in an 
overlap and simulated kymographs of motor trajectories are similar to those found experimentally 

(%[!])• 

We derived analytical and approximate solutions of the continuum steady-state equations and 


compared them to kMC simulation results (fig. S2). Focusing on the symmetric case where the 
boundary conditions a and /3 are the same for both filaments in the overlap, we used both kMC 


simulations and phase plane flows (fig. S3 - S5) to determine the nonequilibrium phases possible 


m 


our model (fig. [^. For low rate of motor switching between filaments in the overlap, we find the 
low density, high density, low density-high density, and Meissner phases previously studied for the 
single-lane case (|5|) • In addition, for high switching rate we find a novel low density-high density-low 
density-high density phase with three domain walls. We determined representative phase diagrams 
for small and large switching rate (fig. . 

We then studied the model in more detail both for the reference parameters and a high- 
switching-rate parameter set (table [STj) . Results of kMC simulations (fig. [S6|) a gree well with 
exact and approximate solutions of the continuum steady-state equations (fig. B . In con¬ 
trast to systems in which motors move on single filaments (|7llI7H22l|2lD, we find that antiparallel 
overlaps with no-flux boundary conditions have total zero current. This leads to a total binding 
constraint that the integral of the total motor density on a single filament must equal po, the motor 


density set by binding/unbinding equilibrium (fig. S8). 


For the experimentally relevant low density-high density coexistence phase, the density profiles 


are approximately piecewise linear (fig. S9). This motivates an analytic approximation to determine 
the density profiles consistent with the total binding constraint and gives analytic expressions for 
the overlap center motor density and the length of the boundary layer in which motors accumulate 
near the overlap ends. We find that as a result of the total binding constraint, the motor density 













Motors on antiparallel MT overlaps 


10 


at the center of the overlap is not determined solely by the motor binding equilibrium, but is also 
controlled by the overlap length, motor speed, and filament switching rate (fig. These same 

parameters control the length of the boundary layer at the overlap ends where motors accumulate 
(fig. [Slop. 

The mitotic spindle contains arrays of overlapping antiparallel MTs to which multiple motors 
and crosslinkers bind (IMp. The surprising differences in motor density prohles between single 
filaments and the antiparallel overlaps we study here are therefore of interest in the study of the 
spindle midzone. For antiparallel overlaps, both the motor density far from the overlap ends and the 
number of motors near the overlap ends can be tuned not just by motor binding kinetics but also by 
motor speed, filament switching rate, and overlap length (hg.[^. The antiparallel filament geometry 
gives biological systems additional handles to control motor density, its spatial distribution, and 
therefore motor function. Motor density can affect recruitment of other proteins and MT dynamics. 
Therefore, this previously undescribed mechanism of regulation of motor density along MTs may 
offer advantages to the control of motor activity. 

For single MTs, length regulation by the kinesin-8 motor Kip3 depends on a functional form 
of motor density that is independent of MT length (|71 [T6]) . The physics we describe here means 
that the motor density and extent of motor accumulation at the overlap ends depends on overlap 
length. Therefore, length regulation of antiparallel MT overlaps must occur differently. In future 
work, it will be of interest to understand how the motor density profiles on fixed-length overlaps 
can be used to understand the length regulation of dynamic overlaps. 
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Figure Legends 

Figure 

Model and results overview. (A) Schematic of the model of motor motion on an antiparallel 
microtubule overlap. Two filaments (green and blue) are modeled as ID lattices with their plus 
ends oppositely oriented. The filaments are labled R (L) if the plus end is pointing to the right 
(left). Motors (red) bind to empty lattice sites with rate konC and unbind with rate kos. Bound 
motors step toward the MT plus end with rate v (if the adjacent site toward the MT plus end is 
empty) or switch to the other MT with rate s (if the corresponding site on the adjacent MT is 
empty). At MT minus ends, motors are inserted at rate av. At MT plus ends, motors are removed 
at rate /3u. (B) Simulated experimental images made from our kMC model. Green, motor density. 
Red, overlap region. Scale bar, 5 pm. Simulations used to generate these images used the reference 
parameter set (table and the indicated bulk motor concentrations. (C) Simulated kymograph 
made from our kMC model with motor spatial position on the horizontal axis and time increasing 
down. Horizontal scale bar, 10 ;um. Vertical scale bar, 5s. The simulations used to generate the 
kymograph used the reference parameter set and 0.5 nM bulk motor concentration. 

Figure 

Nonequilibrium phases. Left: motor density Pr{x) on MT with rightward-moving motors. Right: 
total motor density Pr{x)+pl{x) on both MTs in the overlap. The parameters used are the reference 
parameter set of table with a bulk motor concentration of 200 nM, except for the switching rate 
which is 0.5 s“^ for the LHLH curve and 0.1 s“^ for all other curves and the boundary conditions 
as noted below. LHLH: low density-high density-low density-high density coexistence {a = 0.3, 
j3 = 0.9, orange). Hg: high-density end-accumulating phase (a = 0.95, (3 = 0.99, light blue). He: 
high-density center-accumulating phase (a = 0.6, f3 = 0.95, dark blue). LHc: low density-high 
density center-accumulating coexistence (a = 0.3, f3 = 0.95, green). M: Meissner phase a = 0.6, 
j3 = 0.4, yellow). Lg: low-density end-accumulating phase (a = 0.3, {3 = 0.4, black). LHg: low 
density-high density end-accumulating coexistence (a = 0.1, f3 = 0.95, magenta). Lc: low-density 
center-accumulating phase {a = 0.05, /3 = 0.1, red). 

Figure 

Phase diagrams. Left: low switching rate (0.1 s“^); right: high switching rate (0.5 s“^). Solid 
lines indicate phase boundaries, and dotted lines boundaries between center-accumulating (v) and 
end-accumulating (e) overlap motor density profiles. The bulk motor concentration is 200 nM, the 
motor speed is 5 pm s“^, and other parameters are the reference values (table [^. 

Figure 

Motor density profiles for the reference parameter set. Left: motor density Pr{x) on MT with 
rightward-moving motors. Right: total motor density Pr{x) + Pl{x) on both MTs in the overlap. 
Parameters are the reference parameter set of table with the bulk motor concentrations indicated 
in the legend. 

Figure 

Motor density at the center of the overlap. Left: variation with bulk motor concentration; right: 
variation with motor speed. Points indicate simulation results and dashed lines theoretical pre- 
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diction from Eqn. S31 
parameters. 


The red line shows the value of the Langmuir density pQ for the same 


Figure 

Dependence of density profiles on motor speed. Left: motor density Pr{x) on MT with rightward- 
moving motors. Right: total motor density Pr{x) + Pl{x) on both MTs in the overlap. Varying 
motor speed can significantly alter the motor density at the center of the overlap. These simulations 
use the large-S parameter set of table with bulk motor concentration c = 200 nM and the motor 
speeds indicated in the legend. 


Figure 

Control of boundary layer length. Left: boundary layer length as a fraction of total overlap length, 
shown as a function of bulk motor concentration and motor speed. Right: boundary layer length as 
a function of overlap length for the reference parameter set and varying bulk motor concentration. 
The boundary layer length is the distance at each end of the overlap where significant motor 


accumulation occurs and is determined by Eqn. S30 
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Figure 1: Model and results overview. (A) Schematic of the model of motor motion on an antipar¬ 
allel microtubule overlap. Two filaments (green and blue) are modeled as ID lattices with their 
plus ends oppositely oriented. The filaments are labled R (L) if the plus end is pointing to the right 
(left). Motors (red) bind to empty lattice sites with rate konC and unbind with rate kog. Bound 
motors step toward the MT plus end with rate v (if the adjacent site toward the MT plus end is 
empty) or switch to the other MT with rate s (if the corresponding site on the adjacent MT is 
empty). At MT minus ends, motors are inserted at rate av. At MT plus ends, motors are removed 
at rate /3u. (B) Simulated experimental images made from our kMC model. Green, motor density. 
Red, overlap region. Scale bar, 5 pm. Simulations used to generate these images used the reference 
parameter set (table and the indicated bulk motor concentrations. (C) Simulated kymograph 
made from our kMC model with motor spatial position on the horizontal axis and time increasing 
down. Horizontal scale bar, 10 pm. Vertical scale bar, 5s. The simulations used to generate the 
kymograph used the reference parameter set and 0.5 nM bulk motor concentration. 
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Figure 2: Nonequilibrium phases. Left: motor density Pr{x) on MT with rightward-moving motors. 
Right: total motor density Pr{x) + Pl{x) on both MTs in the overlap. The parameters used are 
the reference parameter set of table with a bulk motor concentration of 200 nM, except for the 
switching rate which is 0.5 s“^ for the LHLH curve and 0.1 s“^ for all other curves and the 
boundary conditions as noted below. LHLH: low density-high density-low density-high density 
coexistence (a = 0.3, j3 = 0.9, orange). Hg: high-density end-accumulating phase (a = 0.95, 
j3 = 0.99, light blue). Hg: high-density center-accumulating phase (a = 0.6, (3 = 0.95, dark blue). 
LHg: low density-high density center-accumulating coexistence {a = 0.3, j3 = 0.95, green). M: 
Meissner phase a = 0.6, (3 = 0.4, yellow). Lg: low-density end-accumulating phase (a = 0.3, 
13 = 0.4, black). LHg: low density-high density end-accumulating coexistence {a = 0.1, (3 = 0.95, 
magenta). Lg: low-density center-accumulating phase (a = 0.05, j3 = 0.1, red). 




Figure 3: Phase diagrams. Left: low switching rate (0.1 s“^); right: high switching rate (0.5 s“^). 
Solid lines indicate phase boundaries, and dotted lines boundaries between center-accumulating (v) 
and end-accumulating (e) overlap motor density profiles. The bulk motor concentration is 200 nM, 
the motor speed is 5 pm s“^, and other parameters are the reference values (table [^. 
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Figure 4: Motor density profiles for the reference parameter set. Left: motor density Pr{x) on 
MT with rightward-moving motors. Right: total motor density Pr{x) + Pl{x) on both MTs in the 
overlap. Parameters are the reference parameter set of table with the bulk motor concentrations 
indicated in the legend. 
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Figure 5: Motor density at the center of the overlap. Left: variation with bulk motor concentration; 
right: variation with motor speed. Points indicate simulation resnlts and dashed lines theoretical 
prediction from Eqn. S31 The red line shows the value of the Langmuir density po for the same 
parameters. 
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Figure 6: Dependence of density profiles on motor speed. Left; motor density Pr{x) on MT with 
rightward-moving motors. Right: total motor density Pr{x) + Pl{x) on both MTs in the overlap. 
Varying motor speed can significantly alter the motor density at the center of the overlap. These 
simulations use the large-5 parameter set of table SI with bulk motor concentration c = 200 nM 
and the motor speeds indicated in the legend. 
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Figure 7: Control of boundary layer length. Left: boundary layer length as a fraction of total 
overlap length, shown as a function of bulk motor concentration and motor speed. Right: boundary 
layer length as a function of overlap length for the reference parameter set and varying bulk motor 
concentration. The boundary layer length is the distance at each end of the overlap where significant 


motor accumulation occurs and is determined by Eqn. S30 




































